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The nonlinear photon-photon interaction mediated by a single two-level atom is studied theo- 
retically based on a one-dimensional model of the field-atom interaction. This model allows us to 
determine the effects of an atomic nonlinearity on the spatiotemporal coherence of a two photon 
state. Specifically, the complete two photon output wave function can be obtained for any two pho- 
■ ton input wave function. It is shown that the quantum interference between the components of the 

f"^ | output state associated with different interaction processes causes bunching and anti-bunching in the 

, two photon statistics. This theory may be useful for various applications in photon manipulation, 

£S| ' e.g. quantum information processing using photonic qubits, quantum nondemolition measurements, 

, t ' and the generation of entangled photons. 

^ . PACS numbers: 42.50.Ct, 42.65.-k, 32.80.-t 

(N ■ 

(NJ ' I. INTRODUCTION 

> : 

| The nonlinearity of atomic objects, e.g. two-level atoms and quantum dots, can be sensitive to individual photons. 
I . This kind of nonlinearity may be useful for the study of photon manipulations such as quantum information processing 
t-H ' IH Q i quantum nondemolition measurements [BJ , and the generation of entangled photons [(| . Realizations of 
single atom nonlinearities have been studied extensively in the field of cavity quantum electrodynamics M, la, 111 • The 
sensitivity of this atom-cavity system to individual photons has been demonstrated by Turchette et al. Q] • Recently, 
we proposed a modification of this setup that enhances the nonlinearity by avoiding all losses in a one-sided cavity 
geometry Every input photon will then be found in the output. The nonlinear response of a single atom to 
i an input of two photons, e.g. from single photon sources |Toj . can then be observed in the correlations between 
^Jj i the two output photons. In order to apply the correlations due to the nonlinear photon-photon interaction, e.g. in 
quantum information processing, it is important to understand the precise spatiotemporal coherence of the input and 
output photons. This cannot be fully achieved by theories that eliminate the quantum state of the field outside the 
atom-cavity system pj yj, |9j . We therefore propose a theory that includes the propa gat ion to and from the system in 
the quantum state, based on a one-dimensional model of the field-atom interaction pd| . In this paper, we apply this 
one-dimensional model of the field-atom interaction to the case of two photon input wave packets. 

The rest of this paper is organized as follows. In sectional we give a theoretical model of the light-atom interaction 
in one-dimensional free space. In section ITTT1 we discuss the experimental realization. In section llVl we derive the 
general solution of the Schrodinger equation for the one photon case. In section |Vj we apply the result of the one 
photon case to derive the general solution for two photons. In section IVT1 the analysis of spatiotemporal coherence in 
the outgoing wave packet is performed as an example for the applications of our theory. In section fVlII the effects of 
the nonlinear interaction on the second order correlations is discussed. In section lYUH it is shown that the bunching 
and anti-bunching effects in the two photon statistics can be understood as quantum interference effects of different 
output components. In section llXl we conclude with a summary of our discussions. 



II. THEORETICAL MODEL OF THE LIGHT-ATOM INTERACTION IN ONE-DIMENSIONAL FREE 

SPACE 

To investigate the change of the spatiotemporal quantum coherence originating from the nonlinear interaction of 
two photons mediated by a single two-level system, we need a model for the spatiotemporal propagation to and from 
the atom. A possible model has been studied in the analysis of spatiotemporal quantum coherence for the case of 
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spontaneous emission from a single excited atom |xT|. This model is illustrated in fig. ^ The r-axis represents the 
single spatial coordinate of the one-dimensional light field. A single two-level atom is coupled locally with the light 
field at the position r = 0. The negative region r < and the positive region r > correspond to the incoming 
field and the outgoing field, respectively. This means that the light field can only propagate in the positive direction, 
approaching the atom at r < 0, and moving away from it at r > 0. The dispersion relation describing the field 
dynamics is then given by the wave number multiplied by the speed of light, u) — ck. The Hamiltonian composed of 

r (space) 

r > 




single two-level atom 



r < 

FIG. 1: Illustration of the theoretical model. The r-axis represents the single spatial coordinate of the field. A single two-level 
atom is placed at the position r — 0. E and G represent the ground state and the excited state of the atom, r > corresponds 
to the output field and r < corresponds to the input field. 



the uninterrupted propagation and the interaction between the atom and the one-dimensional field can be written as 

/oo j -p 

dk ih\ — (bUk)&- -olS(fc)] , (1) 
-oo V 7T V J 



where b(k) is the photon annihilation operator, and tr_ is the annihilation operator of the atomic excitation. \J cT/tt 
is the coupling constant between the atom and the light field. T is the dipole relaxation rate. As will be seen later, this 
rate defines the only relevant time scale of our model. Note that the Hamiltonian has been formulated in a rotating 
frame defined by the transition frequency ujq. Likewise, the wave vector k is defined in the rotating frame, i.e. k is 
defined relative to the resonant wave vector uiq / c. 



III. EXPERIMENTAL REALIZATION 



The situation described by the theoretical model in section [H] can be realized experimentally by using a one-sided 
cavity as illustrated in fig. El The left mirror of the cavity has a transmittance much higher than the right mirror, 
which has nearly 100 % reflectance. The negative region on the space axis in the model shown in fig. ^ corresponds 
to the input in fig. Eland the positive region corresponds to the output in fig. El 

In terms of the conventional cavity quantum electrodynamics parameters, this regime is characterized by k > g, 
where re is the cavity damping rate through the left mirror and g is the dipole coupling between the atom and the 
cavity mode. Therefore the method of adiabatic elimination can be applied to the time evolution of the cavity field Q. 
This means that since the cavity damping rate re is much faster than the dipole coupling g, the interaction between 
the atom and the outside field mediated by the cavity field can be expressed by an effective dipole relaxation rate 
r = g 2 j k. It should be noted that, in this case, the Hamiltonian (1) represents an approximation valid only within the 
finite cavity bandwidth of 2re. Effectively, the theoretical model can then be used to correctly describe the atom-cavity 
dynamics at timescales larger than 1/re. This approximation is sufficient for the description of atomic absorption and 
emission processes as long as l/r> 1/re. Features of the cavity response which become obvious only at a time scale 
of about are neglected in this paper. 
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The dipole relaxation rate T describes the dipole damping caused by emissions through the left mirror of the cavity, 
and the corresponding rate of spontaneous emission through the cavity is equal to 2T 0011 In °ur case, we assume 
that the rate of spontaneous emission into the non-cavity modes 711 is negligible (711 <C 2F). Nearly all emissions from 
the atom can then be confined to the cavity and 2r is the total spontaneous emission rate of the excited atom in the 
cavity. In present cavity designs, this can be realized by covering a large solid angle of the atomic emission with the 
cavity mirrors and exploiting the enhancement of spontaneous emission by the cavity. For example, in the case of 
Turchette et al.'s experiment @, the cavity parameters indicate that about 70% of the sponaneous emission from the 
atom is emitted through the single cavity mode. Another promising method of achieving a one-dimensional emission 
and absorption of the atom is the use of semiconductor microstructures, as reported e.g. in 12] . In any case, our 
model should apply to any cavity design with k 3> g S> 711 . 



T T ~ 




FIG. 2: Schematic representation of cavity geometry. T and T are the transmittances of the mirrors. E and G represent the 
ground state and the excited state of the single two-level atom. The arrows to the left of the cavity represent the free space 
input and output fields, respectively. 



IV. ONE PHOTON PROCESSES 



In this section, we treat the interaction of one photon with the atomic system as a preparation for the analysis 
of two photon processes. We can expand the quantum state of the single photon in the basis of the wave number 
eigenstates \k) and the excited state |E) of the two-level atom. The quantum state for the one photon process can 
then be written as 



|*(t))=*(E;t)|E)+ / dk *(fe;t)|fe). (2) 

-' CO 

In this basis, The Hamiltonian given by Q can be expressed as 

-Blphoton = frck + fli n t (3) 

/CO / -p poo 

dkk\k)(k\ and H iat =ihJ— / dk (\k)(E\ - \E)(k\) . 

"GO V 7T J QQ 

The equations for the time evolution of the probability amplitudes ^(E; t) and ^(k; t) can then be obtained from the 
Schrodinger equation ih d/dt \^{t)) = using © and ©, 

^*(E;i) = -J— / dk*(k;t) (4) 

at v 7r j_ oa 

4 = -ikc^(k;t) + \ — *(E;t). (5) 

dt V 7T 

The time evolution t) can be obtained by integrating 

*(k;t) = e- lkc{t " u) ^(k;ti) + J— / dt e -*M*-* >¥(E;t), (6) 
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where ti is the initial time of the time evolution. In order to describe the time evolution in real space, the result of 
the integration © can be Fourier transformed using 

1 f°° 

^{r;t) = ^= dk e tkr ^{k,t). (7) 



/2tt 

The real space representation of the time evolution then reads 

,'*(r - c(t - ti);ti) for r < or dt - U) < r 

ty(r-t) — I i — (8) 
y,/ \*(r-c(t-^;*i) + \/ir*(E;i-£) for < r < c(t - U). 

The top term corresponds to the single photon amplitude propagating without being absorbed by the atom. On 
the other hand, the bottom term consists of two processes. The first term also corresponds to propagation without 
absorption, while the second term corresponds to the amplitude of a single photon reemitted into the outgoing field 
after absorption by the atom [ijj. The time evolution 4'(E;t) of the excited state amplitude can be obtained by 
integrating (0J using the result for ^(k;t) given in (JSJ, 



*(E;t) = -J— dt dk^{k-t) 

V 7T J t . J_ 00 

' l r dt [ dk (e- ik < t - ti ^{k;t i ) + J— f dt e- lfcc < t -*')*(E;i') ] (9) 



Using the Fourier transform to obtain the real space representation of \t(fc;t), the result reads 

*(E; t) = e - r( *- tl) *(E; U) - VzFc f dt e^^tf (-c(i - U);U). (10) 



The first term corresponds to emission from the excited atom and the second term corresponds to the excitation of 
the atom by the incoming light. 

Since we are interested in the response of the ground state atom to a one photon input, the initial conditions of 
^(Ejii) and $>(r;ti) can be defined as follows, 

V{r>0;U) = 0, (11) 
*CE;t<) = 0. (12) 

Condition (|llfl represents the assumption that there is no light at r > 0. Condition (|12f> represents the assumption 
that the atom is initially in the ground state. With the above conditions, the time evolution ^(r;t) given by iJSJ and 
(fTUf) becomes 

{^{r-c(t-ti)\ti) forr<0 
*(r - c(t - U); U) - 2r f t ; r/c dt e -r(t-r/c-t') $ ( _ c(t ' _ u) . f .) for < r < c(t - U) (13) 
for r > c(t - ti). 

To investigate the outgoing amplitude \&(r > 0;t) for an arbitrary incoming amplitude under the conditions l|llf) 
and iJUJl, it is convenient to represent \&(r > 0;t) by using the matrix element of the time evolution operator. In 
order to derive this matrix element, we expand the time evolution ^(i)) of the quantum state given by (|11|) and 112|) 
as follows, 

!*(*)) = u(t-u)Mti)) 

drdr\r)(r\fi(t-ti)\r){r'\*(ti)) 



+ |E) / dr (E|#(f-t<)|r )(r |tffo)>, (14) 



where |r) = dk e lkr \k) is the eigenstate of a photon at the position r. U(t — ti) = e i H iph°ton (* U) i s the 

time evolution operator. ^(r;i) can then be expressed as follows, 

*(r;i) - (r|*(t)) 

dr Ui photon (r, r ; t - U)^(r ; U), (15) 



where u lphoton (r;r ;t-ti) = {r\U (t - U)\r ). 
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uiphoton(f; r ;t — ti) is the matrix element of the time evolution operator U(t — ti). This matrix element is the 
transition probability amplitude from the state \r ) at time ti to the state |r) at time t. Ui p hoton(?", f ]t — ti) can be 
obtained by comparing the results (|13|l and (|15() . It is given by 



ui p hoton(r; r ;t — ti) 
with 

and 



u P rop(r; r ; t - U) + u abs (r; r ;t - U) 
u P rop(r; r ;t-t i ) = S(r - c(t - U) - r ] 

'_2T e -E(c(t-i i )+r' 





u a b s (r;r ;t- U) = 



for < r < c(t - U 



and r < 



for r > c(t — t^ + r or r > 0. 



(16) 



This transition amplitude can be interpreted as follows. At the initial time t iy the photon starts to propagate at the 
position r . For r < 0, it propagates at constant velocity c. For r > 0, it has already passed the atom and it is in 
a superposition of two amplitudes corresponding to the uninterrupted propagation of the photon and to absorption 
and reemission by the atom, respectively. u prop corresponds to the uninterrupted propagation and u a b s corresponds 
to the reemission from the atom at an emission rate of T. An illustration of the one photon time evolution is shown 
in fig. |21 
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FIG. 3: Schematic representation of the time evolution of incoming and outgoing amplitude. The horizontal axes in (a) and 
(b) represent the space coordinate r. The vertical axes represent the probability amplitude of the spatial one-dimensional field, 
(a) illustrates the amplitude before the arrival of the photon at the atom and (b) illustrates the amplitude after the arrival of 
the photon at the atom. In (a), the photon, having a delta like amplitude, is propagating in the incoming field. The amplitude 
at time t has its peak at r = c(t — ti) + r . In (b), the outgoing amplitude of the photon is in a superposition of two amplitudes 
associated with the uninterrupted propagation of the photon and with reemission after absorption by the atom, respectively. 
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V. TWO PHOTON PROCESSES 



In this section, we treat the interaction of two photons with a single two-level atom. The method we use here to 
describe the two photon quantum state is to first distinguish the particles and then to introduce the correct symmetry 
of the wavefunction for indistinguishable bosons. Note that this is a standard textbook approach to problems in multi 
particle quantum mechanics [13. If we consider the two photons to be separate physical systems, we can describe 
their quantum state in the product space of their single photon Hilbert spaces. A general state |^(t)) for two photon 
processes can then be expanded in the product basis of the eigenstates |fc) and |E) as 



-00 
00 



|*(t)> = / dkx tf(fci;E 2 ;t)|fci}® |E 2 ) + / dk 2 ¥(Ei; fc 2 ;*)|Ei) ® \k 2 ) 

) J — OO 

dfadk 2 *(Ai;fc2;t)|Ai) ® \k 2 ) + *(Ei;E 2 ;t)|E 1 ) ® |E 2 >, (17) 



where the indices 1 and 2 of the eigenstates \k) and |E) distinguish the two photons. Since photons are indistinguishable 
bosons, the actual quantum state has to be a state of positive symmetry. This bosonic property of photons requires 
that the quantum states fulfill the following conditions, 

*(fc;E 2 ;t) = *(Ei;fc;t) 

*Oi;fc 2 ;t) = V(k 2 ;ki;t). (18) 

In order to make it easier to formulate the matrix element of the time evolution for two photon processes, we will 
first divide the Hamiltonian into a linear term and a nonlinear term. The linear term describes the dynamics of two 
photons that are absorbed and emitted independently. The Hamiltonian of this interaction free dynamics is then 
given by the sum of the single photon Hamiltonians. Note that this linear Hamiltonian includes transitions to the 
state |Ei,E 2 ), where both photons are absorbed by the atom. Since this transition is impossible in a two-level atom, 
there will be an interaction between the two photons. This interaction can be described by a nonlinear term that 
suppresses the transitions to the state |Ei, E 2 ). In the following, we first introduce the linear term and then solve the 
Schrodinger equation for the linear system using the results of the one photon case. Next, the nonlinear term is added 
to the linear term and the matrix element of the time evolution for two photon processes is derived by a comparison 
between the linear equations and the nonlinear equations. 

The linear component of the two photon Hamiltonian can be expressed using the single photon Hamiltonian (j^Jl, 

rrlin _ rV(l) „ f(2) f(l) rr(2) 

n ~ lphoton lphoton lphoton "^lphoton' 

/OO 
dk i \k i )(k l \ + |Ej)(Ej| 
-OO 

/OO I xT pOC 

dk i k i \k i ){k i \ and ft®=ifiJ— / dki flfciXE^-lEiXkl). (19) 
-00 » ^* J —00 

This Hamiltonian describes the interaction free evolution of the two photon state. Therefore the time evolution of 
the two photon probability amplitudes ^ (k±,k 2 ;t), ^''"(Ei, k 2 \ t) and ^ lln (fci, E 2 ; t) given by this Hamiltonian 
corresponds to the results of the one photon case given in 1 1151) . It reads 

^* Un (fci,fc 2 ;t) = -* C (fc 1 +fc 2 )v|/ lin (fc 1 ,fc 2 ;i) + ,/^(vI/ li "(fc 1 ,E 2 ;t) + * lin (E 1 ,fc 2 ;t)) (20) 
at V 7r 

-* lin (E 1 ,fc 2 ;t) = -ick 2 y Un (V u k 2 ;t)-J— dk x * lin (fci, fc 2 ; i) + W — * lin (Ei, E 2 ; i) (21) 

dt V 7T J_ 00 V 7T 

^(k^E.A) = -ick^ih^t) - J — dk 2 ^(k u k 2 ;t) + J— ^"(Ei.E,;*). (22) 

dt V 7T J^^ V 7T 



Since the time evolution described by H hn corresponds to the time evolution of the single photon case, the integration 
can be performed according to the procedure in eqs. l|tjllU|) . and the matrix element U2p hoton (ri, r 2 ; r l5 r 2 ; t — U) of 
the time evolution operator can be expressed as a product of the individual single photon matrix elements given by 

u 2 photon( r i^ r 2;r 1 ,r 2 ;t-t i ) = ui p hoton(ri; r x ; t - U) ■ ui p h oton (r 2 ; r 2 ; t - U). (23) 
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The photon-photon interaction can now be included by adding the nonlinear term that suppresses transitions to the 
two photon absorption state |E!,E 2 ), 

tt Tjlin _l_ A fVNonlin 

-H2photon — tl + IXtl , 

where AtfNoniin = _ + |Ei)(Ei| £<2)j (24) 

With this addition, the matrix elements of the two photon Hamiltonian (24) are identical with the corresponding 
two photon matrix elements of the original Hamiltonian (1). The time evolution of the wavefunction is modified only 
slightly by the addition of the nonlinear term. It now reads 

^(k u k 2 ;t) = -ic(k 1 + k 2 )*{k 1 ,k 2 ;t) + J— (*(/Ei,E 2 ;t) + *(Ei,fc 2 ;*)) (25) 
at \ tt 



-^(E^k^t) = -icfc2*(E ls fe2;t)- J— / dh^ikuhit) (26) 
at V tt J_ OQ 

4-(fci,E 2 ;f) = -icfci#(fei,E 2 ;t)-\/— / dk 2 ^(h,k 2 ;t). (27) 



The comparison between <|21l22ll and (|26I27() shows that the integration of 126127ft is the same as the integration 
of I|21I22|I . except that 4 , (Ei,E 2 ;i) is zero in the integration of ^(Ei, k 2 ;t) and ^(fci, E 2 ; t). This means that the 
matrix element u 2p hoton (ri , r 2 ; rj , r 2 ; i — U) of the time evolution described by -ff 2p hoton can be obtained by identifying 
the components in the interaction free propagation of the two photons from r 1 , r 2 to r\ , r 2 described by the matrix 
element u lln (ri, r 2 ; r l , r 2 ; t — ti) and removing the ones due to two photon absorption. 

For the case that photon 1 and photon 2 start at r 2 < r 1 < 0, the time evolution of this probability amplitude can 
be obtained by distinguishing the following three time regions: 

Time region I: t < U — —, the two photons are independently propagating in the incoming field (See fig-E](a)). The 
matrix element associated with the uninterrupted propagation uJjphotonC 7- !' r 2 ; r 1 ,r 2 ;t — ti) can be expressed as 

u l P hoton( r ii r 2]r 1 , r 2 ;t - U) = u prop (ri; r x ; t - U) ■ u pr0 p(r 2 ; r 2 ; t - U) . (28) 

Time region II: t i L < t < ti — — , photon 1 has already arrived at the atom and photon 2 has not arrived yet. In 

this situation, photon 1 is interacting with the atom and photon 2 is propagating in the incoming field (See fig. 01 
(b)). The matrix element associated with one photon absorption u 2 p hoton (ri, r 2 ; r lf r 2 ; t — ti) can be expressed 
as 

u 2photon( r i> r 2;?'i,r' 2 ;i- U) = uiphoton^i;^;*- ti) ■ u prop (r 2 ; r 2 ; t - U) . (29) 

Time region III: ti < t includes a new situation not treated in the one photon processes, because both pho- 
tons have now reached the atom. In this time region, the time evolution given by -ff 2p hoton is different from 
the one given by H hn because there cannot be any contributions of the two photon absorption amplitude 
v I'(Ei,E 2 ;t — ti) in the time evolution of the outgoing amplitude. This means that the matrix element corre- 
sponding to two photon absorption u a b s (r x ; r x ; t — ti) ■ u a b s {r 2 ; r 2 ; t — ti ) will not be included in the time evolution 
of u 2p hoton(?'i, f 2 ] Ti, r 2 ; t — ti) if the position n of photon 1 is less than c(t — U) + r 2 , since, in this case, photon 
2 arrives at the atom before photon 1 is reemitted. The matrix element u 2p I hoton (ri , r 2 ; r 1 , r 2 ; t — ti) can then 
be expressed as 

in , ' ' , ,-. I u^hotoniri;^; t - U) ■ u lpho t on (r2;r' 2 ;t - U) for n > c(t - U) + r' 2 

[uiphotoniri;^;* - tj) • u prop [r 2 ;r 2 ;t - U) for r x < c(t - U) + r 2 . 

Note that the dependence of the transition amplitudes for photon 2 on the output coordinate of photon 1 makes 
it impossible to separate the dynamics of photon 1 and photon 2. Therefore it is not possible to illustrate this 
time region using single photon amplitudes. 
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FIG. 4: Schematic representation of the time evolution of the incoming and the outgoing amplitude in time region I and II. 
The horizontal axes in (a) and (b) represent the space coordinate r. The vertical axes represent the probability amplitude of 
the spatial one-dimensional field, (a) corresponds to time region I, (b) corresponds to time region II. In time region I, photon 
1 and photon 2, having a delta like amplitude, are propagating in the incoming field. Their amplitudes at time t have their 
peaks at n = e(t — ti) + r 1 and r 2 = c(t — ti) + r 2 , respectively. In time region II, only photon 1 has arrived at the atom, the 
outgoing amplitude of photon 1 is as same as the amplitude for the one photon case. Time region III is not shown, since it 
results in a superposition state that cannot be represented by single photon amplitudes. 



The matrix element with r 1 < r 2 can be obtained from the results for r 2 < r 1 by using the positive symmetry 
between photons propagating in free space. The results of these time regions can be summarized as follows, 

u 2p hoton (n. , r 2 ; r[ , r 2 ; * - U) = u%g hoton (n , r 2 ; r[ , r' 2 ; t - ti) + Au^J 1 " (r i , r 2 ; r [ , r 2 ; t - U ) 
u 2photon { r i, r 2;r' 1 ,r 2l t-t l ) = Uip ho ton {r^r^t - U) ■ Uip hoton (r 2 ; r 2 ; t - U) 

AuSjatri.rajri.r^t-ti) = -^e-?^+^ +2c ^)^^) for < n,r 2 < c(t - U) + Min^, r' 2 ], (31) 

where Min^, r 2 ] is the minimum of r x and r 2 . The matrix element Au^™}" of the nonlinear interaction between the 
two photons originates from the impossibility of two photon absorption at the single two-level atom. The remainder of 
the dynamics is identical to the single photon processes. The total output wave function of two photons propagating 
in the one-dimensional field can then be expressed as 

/oo 
dr 1 dr 2 u 2p h ton (n , r 2 ; r x , r 2 \ t - U) ■ *(r lf r 2 ; U). (32) 
-oo 

The output wave function describes the state of the photons propagating in the far field after the interaction with 
the atom. In general, a two photon wave function propagating in one-dimensional free space obeys the relation 
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^( r i, r 2', t) — V?(ri — ct,r 2 — ct;0). Therefore the function of 3 parameters ^(ri,r 2 ;U) can be expressed by the 
function of 2 parameters ^^(fi — cti,r 2 — cti). In the same way, the two photon wave function ^(r^ r 2 ; tt) in 
the outgoing far field can be expressed as V&out^i — ctf,r 2 — ctf). We can then simplify our result (|3"2"|l by the 
transformation to a moving coordinate system, 

7*1 — ct — X\ 

r 2 — ct = x 2 (33) 
In this coordinate system, the output state at an arbitrary time tj can be expressed as 

*out(^i,a;2) = *(ri,r 2 ;t/) 

/oo 
dx Y dx 2 u(xi, x 2 ; x 1 , x 2 ) ■ ^i a (x 1 ,x 2 ) 
-oo 

with x% = r\ — ctf, x 2 = r 2 — ct f , 

x x = r 1 — cti and x 2 = r 2 — cti, (34) 

where u(xi,x 2 ',x 1i x 2 ) is given by 

u(xi,x 2 ;x 1 ,x 2 ) = u lm (xi,x 2 ;x 1 ,x 2 ) + Au Noilhll (xi,X2;x 1 ,x 2 ), (35) 
where u lln (xi, x 2 ; x v x 2 ) — ui pho ton(a;i; x x ) ■ ui p hoton(x 2 ; x 2 ) (36) 

/ / 2r r ' / 

with ui p hoton(iz;; x ) = 8{x — x ) e~~^ x for x < x (37) 

4p2 , , 

and Au^fn,^;!' il) = e — («i+«m-*i-«*) for Xl , x 2 < Min^i, x 2 \. (38) 

&■ 

Note that the linear component u lln (a:i, x 2 ;x l7 x 2 ) of the matrix element u(xi, x 2 ; x 1} x 2 ) represents the time evolution 
of the two photons without interaction. Therefore it can be expressed by the product of two one photon matrix 
elements. On the other hand, the nonlinear component Au Nonlln (a;i, x 2 ;x 1 ,x 2 ) represents the difference from the time 
evolution of the linear component. The output wave function 5'out(^i, x 2 ) can then be expressed as the superposition 
of the linear term and the nonlinear term by substituting the right hand side of (|35|l into (|34fl . 

*out(zi,Za) = ^ u n t (x 1 ,X 2 )+AvI/Nonlin (xi;:C2) (3g) 

/OO 
dxjdxj, u hn (x 1 ,x 2 ]x 1 ,x 2 ) ■^ in (x 1 ,x 2 ) (40) 
-OO 

/OO 
dx x dx 2 Au Nonlin (a;i, x 2 ; x\, x' 2 ) ■ ^{x^x^. (41) 
-oo 

The formulation above can be used to analyze the outgoing amplitude for any arbitrary incoming two photon state. 
In the next section, these results are applied to the case of rectangular two photon input wave packets. 

VI. GENERAL SOLUTION FOR RECTANGULAR INPUT WAVE PACKETS 

In order to investigate the typical properties of the nonlinear photon-photon interaction, a rectangular two photon 
input wave packet is convenient because the homogeneous probability distribution of the two input photons makes 
it easier to understand the change of the correlations between the two photons due to the interaction. We assume 
that the shape of a single two photon pulse prepared in the incoming far field is a rectangle of L in length. Such a 
rectangular two photon wave packet can be written as follows, 



Vm(xi,X 2 ) = ^in{xi) ■ Vin(x 2 ) 



= \ft ° - X - L (42) 
[0 else. 

Note that, in any practical situation, the flanks of a rectangular pulse will not rise and fall infinitely fast. The discon- 
tinuities of the wavefunction should therefore be interpreted as a continuous change of amplitude that is extremely 
fast on a timescale of 1/r. In particular, the flanks of the pulse should be smooth at a timescale of 1/k due to the 
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limitations of the model regarding the description of the cavity dynamics (see section ITTTf) . However, we assume this 
timescale to be so much shorter than 1/T that its effects can be neglected in the following. 

As shown in J3HJ), the output wave packet can be separated into a linear term and a nonlinear term. The linear 
term in the output of the rectangular input wave packet described by ea. (|42[) can be obtained according to eci. (|4U[) . 
It reads 

/oo 
dx x dx 2 v} m {x ll x 2 ;x 1 ,x 2 ) ■ ^ in (x 1 ,x 2 ) = *out(#i) ■ *out(£2) 
-OO 

with &out(x) = / dx ui p hoton(z; x ) • * in (a; ) 



(43) 



The nonlinear term in the output of the rectangular input wave packet described by eq. (|42|l can be obtained according 
to eq.lglj. It reads 

/oo 
dx[dx 2 Au Nonlin (l£i, X2', x\, X 2 ) ■ ^in(x[,X 2 ) 
-oo 

= ~ e ^ +x ^ ^ e -fMax[o,x 1 ^ a ] _ e -^Ly fa ax < l (i = 1,2), (44) 

where Max[0, x%, x 2 ] is the maximum of 0, x\ and x 2 . The output two photon wave packet can then be written as 

*aut(*i,a: a ) = 9^(xi,x 9 )+A*ZSr aa (xi,X2). (45) 




It should be noted that the shape of the output wave packet is defined by the ratio of the dipole relaxation length 
c/r and the length L of the input wave packet. In frequency representation, this means that the shape of the output 
wave packet is sensitive to whether the frequency spectrum of the input wave packet is narrower than the atomic line 
width 2r or not. A particularly simple case for showing the contribution of the nonlinear term to the output wave 
packet can be obtained in the long pulse limit c/r <C L because, in this limit, the linear term given by ea. (|43|l of the 
output wave packet is almost constant for the region < Xi < L — 2c /T (i = 1, 2), that is, 

$^ t (xi,x 2 ) ^ - for0<a; i <L-2c/T (i = l,2). (46) 

On the other hand, the nonlinear term given by ea. H44|) of the output wave packet can be approximated as 

A^X nn (xi,X2) ^ -~e for < Xi <L-2c/T (i = l,2) (47) 

which only depends on the relative distance \x\ — x 2 \ between the two photons. The exponential decay indicates that 
the nonlinear deviation from the linear term only becomes significant in the vicinity of X \ = x 2 . The output wave 
function can then be written as 

V out {x 1 ,x 2 ) ~ I - l e -Sl*i-«*l for0<x i <i-2c/r(i = l,2). (48) 

The contour plot in fig. 0(a) shows an example of the output wave function x E'out(a;i, x 2 ) in the long pulse limit. 
The probability amplitude increases from black to white shading. In this example, we have chosen an input wave 
packet length L = 20c/r which is 20 times greater than dipole relaxation length c/r. The plateau region away from 
x\ = x 2 in fig. [S] (a) has positive amplitude. On the other hand, in the vicinity of x\ — x 2l a valley of negative 
amplitude cuts across this plateau. The shape of this valley can be seen in fig. El (b). Fig.[5](b) is the cross-section of 
the contour plot at x 2 = 10c/r. It should be noted that the shape of the valley is the same for any cross-section x 2 
within the plateau region. The valley is due to the contribution of the nonlinear term in ea . (|44|l which decreases with 
the distance between x± and x 2 . The plateau is the unchanged characteristic feature of the long rectangular input 
wave packet. Fig[S] shows the typical characteristics in the output for a long two photon input pulse. In the next 
section, these typical characteristics are analyzed in terms of two photon statistics. 
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x\ in units of c/T xi in units of c/T 



FIG. 5: Contour plot of the output two photon wave packet for an input wave packet length of L = 20c/F. The horizontal 
axes xi in both figures represent the space coordinate of photon 1, and the vertical axes X2 represent the space coordinate of 
photon 2. (a) is a contour line plot of ^out^i, £2) and (b) is the cross section of the contour plot at X2 = lOc/F 

VII. ANALYSIS OF TWO PHOTON STATISTICS 

Two photon statistics are often used to characterize the properties of nonclassical light from the viewpoint of photon 

(2) 

counting. For this purpose, the second order correlation G Q ^ t (t,t + r) is defined by the joint probability density of 

detecting one photon at time t + r when the other photon is detected at time t. When the second order correlation 

(2) 2 
G<; ut (i, f+r) is expressed by using the probability density l^outC^ + ct, x)\ , the time t is given by the space coordinate 

x divided by the speed of light, t = —x/c. Also, the probability density |$ ut(a; + ct, x)\ 2 must be rescaled by a factor 
of c 2 , since the probability density per unit time is expressed by c times the probability density per unit length and, 

(2) 

for second order correlations, this factor is squared. Moreover G„ ut does not distinguish between the two photons. 

Therefore, both \ty out (x,x + ct)\ 2 and \^ ou t(x + ct,x)\ 2 contribute to G^ t (t,t + r). Taking all these factors into 
account, the second order correlation function can be written as 

G {2 ^(-x/c,-x/c + t) = c 2 (|*out(a; + cr,x) \ 2 + |* ou t(x, x + ct) | 2 ) = 2c 2 |* out (x + ct,x) \ 2 . (49) 

Note that l^out^ + ct, x)\ 2 and l^out^, x + ct)| 2 always have to be equal to each other because of the bosonic nature 
of the two photons. The contour plot of the second order correlations is shown in fig-H^a-l) and (a-2) for an input 
pulse length of L = 20c/T. The correlation increases from black to white shading. To indicate the contribution 
of the nonlinear term described by eg . (|44|l in the second order correlations, it is convenient to compare this result 
with the second order correlations obtained from the linear component of the output only. Fig. (b-1) and (b-2) 
shows this second order correlation G^ ln (—x/c, —x/c + r) = 2c 2 (^^(a; + ct,x)| 2 of the linear component. The 
comparison of these figures shows that both cases are identical except for their distributions around t — 0. As can be 
seen in fig. (a-2), the nonlinear interaction causes photon bunching around t = and photon anti-bunching around 
t = ±21og2/r. 

In order to compare the second order correlations with other systems in quantum optics, it is useful to normalize 

(2) 

the correlation function G„ ut by the product of the probability densities of single photon detection at times —x/c and 
—x/c + t in the output field. In the long pulse limit, the total probability of finding both photons within the output 
region of < Xi < L (i = 1, 2) is nearly equal to 1. Therefore the single photon detection probability density per 
unit time within this region is 2c/ L. By using this average photon density, we obtain the normalized second order 
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FIG. 6: Distribution of the second order correlations (a) G^ t (— x/c, —x/c + r) and (b) G^ t lln (— x/c + r, —x/c) for an input 
wave packet length of L — 20c/r. The horizontal axes —x/c + t in (a-1) and (b-1) represent the detection time for one photon, 
and the vertical axes —x/c represent the detection time for the other photon. The second order correlations increase from black 
to white shading, (a-2) and (b-2) are the cross-sections of (a-1) and (b-1) at t — 10/r. 



correlation function, 

g™(-x/c,-x/c + T) = G^(-x/c, -x/c + t)/ (2c/ L) 2 

I? 

= y |*out(z + cr,z)| 2 . (50) 
With the approximations for the long pulse limit given by eq.(05J, this correlation reads 

g^ t {-x/c, -x/c + t) ~ - (l-4e- r|r| ) 2 for < {x + ct, x} < L - 2c/T. (51) 

That is, the second order correlation g^J t (—x/c, —x/c + t) in the long pulse limit only depends on the delay time r in 

the output region < {x, x + ct} < L — 2c/T. g^ t around r = is shown in fig. [7] A typical feature of the correlations 
is the double dip feature with zero two-photon coincidence occuring at nonzero time delays of r = ±21og2/r. Such 
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a nonlinear effect is similar to the nonclassical effect which has been reported in [7],ll4j. However, in our case, g^ t 

(2) 

approaches 1/2 beyond the double dip feature. The value g^ ut = 1/2 is a statistical property of the single mode two 
photon input state. Note that the value of g^J t — 1/2 does not indicate that there is a correlation between the two 

(2) 

photons in the pulse. The anti-correlation expressed by g ox ± t = 1/2 simply arises because one cannot detect the same 
photon twice. In the next section, we trace the second order correlation back to the coherent terms of the output 
wave packet. 




-4 -2 2 4 

t in units of 1/r 



FIG. 7: Normalized second order correlation </* u { in the outgoing two photon wave packet. A typical feature of the correlations 
is the double dip feature with zero two-photon coincidence occuring at nonzero time delays ±21og2/F. approaches 1/2 

beyond the double dip feature, since g^ t = 1/2 is a statistical property of the single mode two photon input state. 



VIII. EXPLANATION OF THE SECOND ORDER CORRELATION FUNCTION IN TERMS OF 

QUANTUM INTERFERENCE EFFECTS 

The cause of the features of the second order correlation shown in fig. (a-2) and fig. [7| can be understood 
by analyzing the interaction processes of two photons and the atom as follows. As shown in ea. l|16[) . the matrix 
element of the time evolution Ui p hoton(^; x ) can be expanded in terms of two interaction processes: (I) single photon 
transmission (reflection in fig[2J without absorption; (II) single photon reemission after absorption. In the same way, 
the matrix element of the time evolution u(xi, X2] x 1 , x 2 ) can be expanded in terms of three interaction processes: 
(i) two photon transmission without absorption; (ii) one photon transmission without absorption and one photon 
reemission after absorption; (iii) two photon reemission after absorption, 

u(xi,x 2 ;x 1 ,x 2 ) = \i^\xi 1 X2\x 1i x 2 ) + v^ n \xi,X2\x 1: x 2 ) + vl { -™\xi 1 X2]x 1 ,x 2 )- (52) 

For x 1 < x\ and x 2 < x 2 , the components of u(xi, X2] x^, x 2 ) read 

u( 1 \x 1 ,x 2 ;x 1 ,x 2 ) = u prop (xi; x ± ) ■ u prop (a;2; x 2 ) 

■a. {n \x 1 ,X2]x 1 ,x 2 ) = Uprop^i;^!) • Uabs^;^) + u abs (a;i;a; 1 ) • u piop (x 2 ; x 2 ) 
u (m) (a;i,a:2;x 1 ,a;2) = u abs (xi, x x ) ■ u abs (x 2 ; x 2 ) + Au™ onlm (xi, x 2 ; x 1: x 2 )- (53) 

Therefore, the output wave packet 'I'out^i, X2) can be expanded as 

foutOr 1,^2) = ^l t (xi ) x 2 ) + ^ t (x u x 2 ) + ^(x 1 ,x 2 ) ) 

where V&W to C m ) is given by 

/oo 
dx'.dx'^ to ^(x 1 ,X2,X 1 ,x'2) ■ ttinte'i,^). (54) 
OO 
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For the rectangular input wave packet, the output wave function in the interval < X{ < L (i = 1, 2) is given by 



1 

I 



^ t (x u x 2 ) = I (2e"^-0 _ 2 ) + i (2e-^— J - 2 

a*) - ^(e-i^-^-^e-^-^-lj+A* 1 "" 11 -^!,^). (55) 

Fig. El (a) shows the correlations of ^out* ^(^lj^a) at ^2 = 10c/r for an input pulse length of L = 20c/T. The 
dotted line shows the amplitude corresponding to process (i), where both of two photons are transmitted without 
absorption by the atom. Likewise, the short dashed line shows the amplitude corresponding to process (ii) , where one 
photon is absorbed and then reemitted and another photon is transmitted without absorption. The thin line shows 
the amplitude corresponding to process (hi), where both photons are absorbed and then reemitted. Fig. [S] (b) shows 
the superposition of all the amplitudes. This superposition is equal to the output amplitude VP O ut0&i> x% = lOc/T). 
Therefore fig. El (b) is identical with fig. El (b). Note that the discontinuities of the output wavefunction ^out at = 
and Xi = 20c/r originate from the direct transmission of the rectangular input wavepacket in ^Q Ut and \&" ut . As 
mentioned in section lVll these discontinuities represent changes of the amplitude that are extremely fast on a timescale 
of 1/r, but would be continuous at a much shorter timescale of 1/k. They therefore approximately represent the 
linear response of the atom-cavity system at time resolutions sufficiently longer than the timescale 1/k of the cavity 
dynamics. 

In the long pulse limit L ^> c/T, the amplitudes described by ea. H55|) can be approximated as 



*St(*li*2) 



1 

I 

4 



where A* Nonlin (xi, x 2 ) = — _ e — o'* 1- (56) 

To understand the details of the two time correlations, we now examine the effect of the non-linear contribution 
A^^° n (a;i, x 2 ) in the superposition. The non-linear contribution only depends on the relative distance \x 2 — X\\ 
which corresponds to the difference between the detection times of the two photons. When the relative distance 
\x 2 — x\ \ is much larger than the dipole relaxation length c/T, the non-linear contribution A^E'^J l ° nlln is nearly equal 

to zero and the double absorption amplitude ^out is close to 4/L. On the other hand, when the relative distance 
\x 2 — x\\ is much smaller than the dipole relaxation length c/T, the non-linear contribution A$^ u ° nlln is nearly equal 
to —4/L and the double absorption amplitude ^out is dose to zero. As noted previously, A^™° nlln represents the 
absence of simultaneous double absorption. The time dependence of A4'^ 1 ° nlm therefore describes the saturation 

dynamics of the two-level atom. Since the no-absorption amplitude and the single absorption amplitude 
are independent of the saturation and their amplitudes are 1/L and —4/L respectively, the total amplitude ^out is 
obtained by adding a constant value of — 3/L to ^ut- The resulting total amplitude V&out is then close to 1/L for 
relative distances \x 2 — x\\ ^> c/T and drops to —3/L for relative distances \x 2 — x%\ <C c/T. The total amplitude 
—3/L for relative distances \x 2 — x\\ <§C c/T is associated with the bunching effect in the second order correlations. The 
total amplitude then changes from positive to negative values continuously, depending only on the relative distance 
between the two photons \x 2 — xi\. Precisely speaking, the total amplitude ^out is positive if the relative distance 
is \x 2 — xi\ > 21og(2)c/r and negative if the relative distance is \x 2 — xi\ < 21og(2)c/r (see fig. El (b)). That is, in 
the region with \x 2 — x%\ < 21og(2)c/r, the interaction of the two photons causes a phase flip of tt. This phase flip 
can be understood as evidence on the quantum level for the resonant nonlinearity discussed in Q . The anti-bunching 
at \x 2 — x%\ = 21og(2)c/r originates from ^ ut passing through zero as the sign of the total amplitude changes from 
negative to positive. In this way, both the bunching and the anti-bunching effects in the two photon correlation can be 
explained in terms of interference effects between the quantum coherence contributions from the different interaction 
processes. 
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xi in units of c/Y x\ in units of c/Y 



FIG. 8: Analysis of the interaction processes of two photons and the atom for an input pulse length of L = 20c/F. In (a), 
the dotted line shows the amplitude at X2 = Wc/Y corresponding to process (i), where both photons are transmitted without 
absorption. The short dashed line shows the amplitude corresponding to process (ii), where one photon is absorbed and then 
reemitted while the another photon is transmitted without absorption. The thin line shows the amplitude corresponding to 
process (iii), where both photons are absorbed and then reemitted. The superposition of all the amplitudes is shown in (b). 
Note that (b) is identical with fig. El (b). 



IX. CONCLUSIONS 

We have presented a fully quantum mechanical model of the non-linear interaction of two photons at a two-level 
atom. The experimental realization of such an interaction can be implemented using a one-sided cavity and single 
photon sources. Our theory allows us to determine the effects of an atomic nonlinearity on the spatiotemporal 
coherence of a two photon state. By applying the general results to rectangular input wave packets, we have shown 
that the nonlinear interaction of two photons at the atom causes bunching and anti-bunching effects in the two photon 
output state. Since our model describes the complete spatiotemporal coherence of the output, it is possible to analyze 
the origin of these effects in terms of quantum interference between different absorption and propagation processes. 
This method may therefore provide a useful tool for various applications in the manipulation of individual photons 
such as quantum information processing, quantum nondemolition measurements, and entangled photon sources. 
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